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^ ! Abstract 

o . 

■ The existence of Friedmann limits is systematically investigated for all the hyp- 

. ersurface-homogeneous rotating dust models, presented in previous papers by this 

I author. Limiting transitions that involve a change of the Bianchi type are included. 

Except for stationary models that obviously do not allow it, the Friedmann limit 
expected for a given Bianchi type exists in all cases. Each of the 3 Friedmann models 
qh| has parents in the rotating class; the k = +1 model has just one parent class, the 

^ ' other two each have several parent classes. The type IX class is the one investigated 

bj)' in 1951 by Godel. For each model, the consecutive limits of zero rotation, zero tilt, 

zero shear and spatial isotropy are explicitly calculated. 

; I. Motivation and summary of the method. 

In previous papers^"^ a complete set of all metric forms was derived that can represent 
hypersurface-homogeneous rotating dust models. For each case, the generators of the 
symmetry algebra were found, the Bianchi type determined, and the metric form resulting 
from the Killing equations was explicitly presented. That classification was more detailed 
than the Bianchi classification because all possible orientations of the symmetry orbits in 
the spacetime were allowed, i.e. the orbits could be spacelike, timelike or null. 

In a later paper"^, one of the Bianchi type V models was investigated. Among the 
problems considered there was the question whether the model can reproduce the k = —1 
Friedmann model in the limit of zero rotation, — > 0. Since the coordinates that are 
well-suited to the classification are not suitable at all for considering the limit — 0, this 
limit could be taken only after a coordinate change and reparametrization of the metric. 

In the present paper, the existence of the Friedmann limits is systematically investigated 
for all the other cases found in the classification in Refs. 1-3. The Bianchi type is allowed 
to change in the limiting transition. In all Bianchi type I cases the velocity field is tangent 
to the symmetry orbits, i.e. those models have matter density constant along the flow, and 
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no expanding Friedmann model can be a subcase there. The same is true for the Bianchi 
type II from Ref. 1 and for both the subcases of case 1.1.1.2 in Ref. 2 which are of type 
III. In all the other cases the Friedmann limits that can be expected for a given Bianchi 
type do indeed exist. 

The specialization to the Friedmann metrics is possible in so many cases because there 
is a free parameter in them that determines the tilt of the orbits with respect to the velocity 
field (with various values of the tilt parameter, the orbits may be spacelike, timelike or 
null). Whenever a Friedmann limit exists, the orbits are made orthogonal to the velocity 
field ("untilted") during the limiting transition. 

In order to make this paper readable independently of the other ones, the basic facts 
are briefly recalled here. More details can be found in Ref. 1. 

The velocity field of a rotating dust, m", defines 3 scalar functions T{x),r]{x) and ^{x) 
such that: 

Ua = T,a+r]^,a- (1-1) 

These functions (whose existence follows from the equations of motion via the Darboux 
theorem^) are determined up to the transformations: 

r = r'-5(r,r?'), ^^F{e,v'): V ^ G{i\r^'l (1.2) 

where the functions F and G obey: 

F,^/ G,n' -F,jj' G,^' = 1, (1.3) 

and then S is determined by: 

S,^, = GF,^, -rj', S,,^ = GF,^, . (1.4) 

(eq. (1.3) is the integrability condition of (1.4)). 

Then, the continuity equation, {nu°')-a = 0, where n is the number density of the dust 
particles, implies that there exists one more function ({x) such that: 

V^nu^^e''^^'^,pri„C,s, (1.5) 

where g is the determinant of the metric tensor and £°'l^'^^ is the Levi-Civita symbol. The 
function ( is determined up to the transformations: 

C^C + Tie,v'), (1-6) 

where T is an arbitrary function. 
The following relations hold then: 



9(r,7/,e,C) 



(1.7) 



d{x^, x^,x'^, x^) 

This shows that {T,^,ri,(} can be chosen as coordinates, with r being the time coordi- 
nate. They are called the Plebaiiski coordinates. Denoting {T,^,ri,(} = {x^ , x^ , x"^ , x^} — 
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{t,x,y, z}, we obtain for the velocity field u", the metric tensor Qa/s, the rotation tensor 
U!ai3 and the rotation vector w" in these coordinates: 

S'oo = 1, 901 = y, 902 = 903 = 0, g = det{gai3) = -n~'^ , 

W" = a;c,/3 = -OJ/Sa = ^5^a5^/3- (1-8) 

It is the last property that makes the limitng transition a; — > impossible without a 
coordinate transformation and reparametrization. 

In these coordinates, if any Killing field is allowed by the metric it must be of the form: 

A;" = (C + - )5\ + (P,y 5^1 - 0,, 5% + X5^s, (1.9) 

where C is an arbitrary constant and 4>{x, y) and A(a;, y) are arbitrary functions. If (p^a 7^ 
(i.e. (f) is not constant on an open set), then the coordinates can be adapted to k"" within 
the Plebanski class (by eqs. (1.2) - (1.4) and (1.6)) so that: 

A;" = (5°i. (1.10) 

The metric then becomes independent of x, and the coordinates preserving (1.10) are 
determined up to the transformations: 

t' = t- JyH,ydy + A, x' = x + H{y), y' = y, z' = z + T{y), (1.11) 

where A is an arbitrary constant and H, T arc arbitrary functions. 

If (p^a = 0, then the form of the Killing field /c" = CS°'o + A5"3, is invariant under (1.2) 
- (1.6) and the Plebanski coordinates cannot be adapted to k°'. The property = const 
is equivalent to the following invariant relation: 

k"" = Cii" + (A/n)«;", (1.12) 

i.e. /c" is then spanned on the velocity field and the rotation field. 

If three Killing fields exist, then each of them can either be of the special type (1.12) or 
of the general type (1.9). One of the general-type Killing fields can always be transformed 
to the form (1.10) by (1.2) - (1.6). 

This observation gives rise to a complete classification of all hypersurface-homogcneous 
spacetimes that are possible for a rotating dust. When all 3 Killing fields are of the special 
type (1.12), the symmetry orbits are two-dimensional, and this case is not considered. 
When two Killing fields are of the special type, while the third one is general, there exist 
two classes of metrics (Bianchi types I and II) that were derived in Ref. 1. When one 
Killing field is of the special type, while the two others arc general, all Bianchi types 
except VIII and IX appear (Ref. 2). When all 3 Killing fields are of the general type, all 
the Bianchi types appear, some of them hidden as limits of more general types (Ref. 3). 
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The multitude of cases is a consequence of the many possible alignments or misalignments 
among the 3 Killing fields and the velocity and rotation fields. 

When the Bianchi classification is introduced, the generators of symmetry are scaled 
to standard forms such that all nonzero structure constants (except the free parameters 
in types VI^ and Vllft) become equal either to +1 or to —1. In general, though, they 
arc arbitrary constants, and in the general form each of those constants can be allowed to 
become zero. In this way, the more special Bianchi types can be obtained from the more 
general ones by going to the zero limit with some of the structure constants. The resulting 
hierarchy of Bianchi types is well-known, and is shown in Fig. 1 (adapted from Ref. 5) 
for easy reference. The specializations that are possible can be instantly guessed from the 
values of the a,ni,n2 and parameters for the different Bianchi types. Type III cannot 
be specialized to IV or V because, with the arbitrary values of the parameters n2 and n^, 
the parameter a is determined by a = ^f—n^jiz. 

Another well-known result^ is the placement of different Robertson- Walker geometries 
within the Bianchi classes. This is also recalled for easy reference. Since we are considering 
only dust models, we will call these geometries the Priedmann models and Priedmann limits 
of the rotating models. 

The /c = model is a subcase of the Bianchi types I and VIIq (the two Bianchi algebras 
have different bases, but share common orbits). 

The k — —! model is a subcase of the Bianchi types V and VII/j. 

The k = +1 model is a subcase of the Bianchi type IX. 

When considering each case of the classification from Rcfs. 1-3, one has to recognize 
from Fig. 1 which of the four types {I, V, VIIq, VII/^, IX} could possibly be contained 
in it as a subcase and then the appropriate specialization of the arbitrary constants and 
functions in the model has to be considered. This procedure will be presented in more 
detail in sec. II, later it will be applied without detailed explanations. 

It will turn out that only the stationary models have no Priedmann limit. In every 
nonstationary case, the Priedmann limit indicated by the Bianchi type indeed exists. Note 
that the limits are found for the metrics, without taking into account the Einstein equa- 
tions. This is why a nonstationary type II metric exists in the collection, and is found to 
admit the k = Priedmann limit, even though it is known ^'^ that spatially homogeneous 
type II dust solutions must have zero rotation, see sec. III. 

Now we shall systematically go over all the cases presented in Refs. 2 and 3. The two 
cases from Ref. 1 are immediately seen to admit no Priedmann limit: in both of them, the 
velocity field of the dust is spanned on the KiUing fields (see eqs. (7.7) - (7.8) is Ref. 1), 
so the particle number density n will obey n^au"" = 0. Hence, these cases cannot contain 
any expanding Piredmann model because in the latter ri aU'^ 7^ 0. 

Each of the models presented in Rcfs. 2 and 3 that allows a Priedmann limit will be 
first transformed to the Plebanski coordinates (most of them were found in coordinates 
adapted to the KiUing fields that are not in the Plebanski class). Then, each model will be 
transformed by a coordinate transformation and reparametrization of the metric functions 
and constants to such a form in which the limit of zero rotation can be calculated explicitly. 
Then, the Priedmann limits will be calculated by consecutively imposing on the metric the 
conditions of zero rotation, zero tilt, zero shear and spatial isotropy (i.e. constant curvature 
in the 3-spaces orthogonal to the dust fiow). This last condition is not superfiuous, even 
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though dust with zero rotation and zero shear must be a Friedmann model in consequence 
of the Einstein equations^. It is conceivable that no Friedmann limit would exist at all 
in some classes. However, this does not happen, and a spatially isotropic subcase will be 
found to exist in all cases. The corresponding limits of the KiUing fields, where nontrivial, 
will be also calculated and the Bianchi type of the limit determined. 

Since on each of the underlying manifolds five vector fields exist (velocity, rotation and 
the 3 Killing fields), the 5 vectors must be linearly dependent at each point. This linear 
relation allows to identify in each case the parameter that determines the tilt of the velocity 
field with respect to the symmetry orbits - see sec. V. It turns out that this tilt parameter 
is always simply proportional to that defined by King and EUis^. 

II. The cases 1.1.1 of Ref. 2. 

We begin with case 1.1.1.1., which is of Bianchi type III. 

The transformation from the coordinates used in eq. (2.18) of Ref. 2 (that were adapted 
to the Killing fields) to the Plebahski coordinates is given by eq. (2.16) in Ref. 2 (where 
{t',x',Y,Z} are the coordinates of (2.18), and {t,x,y,z} are the Plebanski coordinates). 
The transformed metric is: 

900 = 1, 901 = y, 9o2 = 9o3 = 0, 

912 = hi2 - ^bX3Yh22, 9i3 = C's^is - ^b'^C3Yh23, 

922 = A3^/i22 , 923 = ^CgAa/ias, 933 = {bCsyh33, (2.1) 
where b, C3 and A3 are arbitrary constants, Y and Z are given by: 

Y = -bXst + Aay + bCsz, Z = bX^t + X^y - bC^z, (2.2) 

and all the hij, i,j — 1,2,3 are arbitrary functions of Z. The first fine of eq. (2.1) will 
be the same in all the other metrics transformed to the Plebanski coordinates, so it will 
not be repeated from now on. Since the argument of hij is determined (by the Killing 
equations) only up to a constant factor, we are allowed to rescale it by an arbitrary factor. 
For considering the limit a; — > 0, it will be convenient to assume that the argument of hij 
is: 

T := Z/{bX^) ^t + y/b- {C^/X^)z. (2.3) 

This presupposes that 6 A3 7^ 0, but this condition is included in the definition of case 
1.1.1.1. The limit A3 = can be taken into account after a simple reparametrization, and 
it leads to a stationary solution. The subcase 6 = is degenerate, and it belongs to the 
1.1.2 family. 
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As seen from the last formula in (1.8), the simplest way to let a; — > is to transform y 
as follows: 

y = ^oV, (2.4) 

and then let a;o — > 0, so that the only nonzero component of rotation in the new coordinates 
becomes: 

10^2 = <^oy > 0- (2-5) 

Then, however, the components g'12, g'22 fl'23 of transformed metric would simulta- 
neously go to zero, and the metric would become degenerate {g = 0). To avoid this, /i22 
must be rescaled as follows: 

/l22 = H22/0Jq^- (2.6) 

Then g[2 — 61^0^12 — \b\zYH22/ojQ would become infinite in the limit ujq 0. To avoid 
this, hi2 must be reset so that the infinity is cancelled. Since all hij depend on T", not on 
y, this can be done as follows: 

/112 = H^2/uJo - \{h\3fTH22/uJo'. (2.7) 

The first term in (2.7) contains the ujq in the denominator for greater generality, so that 
g'^2 > H12 7^ 0. Then, to cancel the infinities in g^n, the function hu must be reset 

Wo— »0 

as follows: 

/ill = H^i - \h\X3TfH22luJ0'' - h^Th^2. (2.8) 

The reparametrization (2.4), (2.6) - (2.8) would be sufficient to make the hmit a;o — > 
of the metric (2.1) nondegenerate. However, the hypersurfaces t = const, that become 
orthogonal to the velocity field 11" in the limit cjq ^ 0, would not yet coincide with the 
hypersurfaces of constant matter density. In the Plebahski coordinates, as seen from (1.8), 
the matter density obeys g = —n~^, and so n would depend on (t — C3Z/X3) in the limit 
loq — > 0, i.e. the model would still be tilted. To "untilt" it, we must let C3 — > 0, and this 
requires at least one more rescaling. It will be convenient to redefine C3 as follows: 

C3 = tooc, (2.9) 

so that the untilting occurs simultaneously with cu ^ 0. Then we must rescale h^a'- 

= Hss/uo''. (2.10) 

For greater generality, we will also rescale h23: 

h23 = ^23/c^o', (2.11) 

and then hi^ must be reset as follows: 

hi3 = H,3/u;o - ]f\3TH23/i^Q^. (2.12) 

The transformation (2.4), applied to (2.1) together with all the subsequent reparametriza- 
tions, results in the following metric: 
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500 = 1, 501 = i^oy, 902 = 5-03 = 0, 

511 = TT-^[^'^3^(-2A3y + 2hcz)ujQ + {X^y + hcz){2>X3y - hcz)ujQ^\ 
4A3 

-^{htf + //n - 2hyH^2 + (bXsyfH22, 

9i2 = H12 - bX3^yH22, gi3 = cHi3 - b'^cX3yH23, 
922 = X3^H22, 923 = ^cAgi^aa, 533 = ibcfH33, (2.13) 

where the Hij depend only on t. Here, similarly as in (2.1), the first line will be the same 
for every metric, and so it will not be repeated from now on. 

The metric (2.13) still has nonzero shear. If a Priedmann model is to result from it, 
the shear must be set to zero. The coordinates {t, x, y, z} in (2.13) are now comoving and 
synchronous, so zero shear means that: 

gij = Gijix,y,z)R^{t), (2.14) 

i.e. all the components of the metric must depend on time only through the same factor 
E?{t). This means: 

other Hij{t) = -CijE?{t), (2.15) 
where Cij are unknown constants. With no loss of generality, it may be assumed that: 

C33 = 1. (2.16) 

The metric (2.13) with Hij as in (2.15) - (2.16) will represent a Fricdmann model when 
the hypersurfaces t = const are spaces of constant curvature. In order to calculate this 
curvature, it is convenient to introduce the new constants Du, D12 and D22 by: 

D22^ '■— C22 — C23' -I D12 '■— (C12 — Ci3C23/b)/{X3D22), 

Dn' := Cn - Cn'/b' - L>i2'. (2.17) 

The correct signs for and 1)22^ are guaranteed by the signature of (2.13). Then (2.13) 
may be written as follows: 

ds^ = dt^ - {DuRdx f - R\Di2 - bX3D22y)dx + X3D22dy]^ 

-R\Ci3/b - bX3C23y)dx + AaCsady + bcdz]^ (2.18) 

and the curvature tensor for the spaces t = const may be easily calculated using the 
orthonormal set of differential forms suggested by (2.18). The curvature tensor is: 

R^\, = ^F'G' + F^ R'\3 = it:''23 = —F^G', (2.19) 
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where: 

F:=b/{DuR), G -.^023/022. (2.20) 

The Riemann tensor (2.19) will represent constant curvature when R^^i2 = R^^is- This 
implies 6 = 0, which seems to be a singular limit of (2.18). However, the limit 6^0 may 
be easily incorporated into (2.18) by the following reparametrization: 

Ci3 = D,3b, c = C/b. (2.21) 

After this, the Riemann tensor of the space t — const still has the same form (2.19) - 
(2.20). With 6 = 0, R'^ki = 0, i.e. (2.18) represents then the k = Friedmann model. 
This is the Friedmann limit of the metric (2.1), as expected for Bianchi type III. 

In this case, the coordinates of the Friedmann limit are similar to those usually used 
(they are the nonorthogonal Cartesian coordinates for the flat space t = const). This will 
not be so in most other cases - the coordinate representation of the resulting Friedmann 
limit will be rather exotic, and calculating the Riemann tensor of the subspace t — const 
will be the simplest way to check that it is the Friedmann metric indeed. 

The Killing fields for the metric (2.1) are (see Ref. 2): 

k^i) = S^, kf2) = e^^(<5o" - b5^), kf,^ = C^S^ + X^S^ (2.22) 

As seen from Fig. 1, the algebra of type III can be specialized only to types II and I, 
and so the A; = Friedmann limit is the only one of the three that can be expected 
here. Note that the Killing field k(^2) will have a meaninigful limit cu — > 0, 6 — * only 
if the two limits are tuned so that uio/b > (for example, b = B^/uJq). Then 

1%. := {ujQ/b)k%. > 5^2, which is indeed a Killing filed of (2.18) with 6 = 0. The 

algebra {^(i), /(2), ^"(3)} becomes then Bianchi type I when ujq = 0, as expected. 

The reasoning behind the reparametrizations, and the subsequent calculation of the 
limits of zero rotation, zero tilt, zero shear and constant curvature of the spaces t — 
const, follows the same scheme in all the other cases. Therefore, it will be presented 
in less detail from now on. In some of the cases, the reparametrization that untilts the 
limit a; — s> is a necessary condition for cancelling the infinities introduced by the earlier 
reparametrizations. 

The other two subcases of case 1.1.1 in Ref. 2, i.e. cases 1.1.1.2.1 (eq. (3.16)) and 
1.1.1.2.2 (eq. (3.32)) are immediately seen to allow no Friedmann limit. For both of them, 
the Killing fields are given by (2.22) above with A3 = 0. As seen from (1.8), the Killing field 
fc(3) coincides then with the velocity field of dust, and so both these models are stationary. 

In fact, the last of (2.22) is the linear relation among the five vectors mentioned at the 
end of sec. I, becuase it is equivalent to (1.12). Since with A3 = the velocity becomes one 
of the Killing fields, i.e. becomes tangent to the symmetry orbits, A3 is the tilt parameter. 
More on this - see sec. V. 



III. Cases 1.1.2 of Ref. 2. 

The case 1.1.2.1 is again of Bianchi type III. The transformation back from the coordinates 
of eqs. (4.6) in Ref. 2 (adapted to two Killing fields) to the Plebahski coordinates is given 
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by (4.4) in Ref. 2, with the roles of {x°'} and {x'"'} interchanged. The transformed metric 
is: 

gn = -2{c/a)y - {c/af + r'/in - 2{c/a)\^Yh^^ + {cX^I afhz:^, 
CA3 

912 = hi2 77^23, 9\z = C^Yhi2, - (c/a)A3C3/l33, 

aY 

922 = h22/Y'^, = C^h2z/Y, 5(33 = (73^/133, Y := ay + c, (3.1) 

where a, c and A3 are arbitrary constants and hij are arbitrary functions of the variable: 

T:^t- C^z/X^. (3.2) 

The reparametrization that will allow setting the rotation and tilt to zero is: 

(y,C3) ^u;o(y,D), 

hu = Hii + (Xs/afH^s/uo'^, = H12/UJ0 + {X^/a)H2z/uJo^ , 

hi3 = Hi3 + {Xs/ajHss/uo^, (/122, /i23, hs) = (^22, -f^23, H33)/ujQ^. (3.3) 
The reparametrized metric (without the limit a;o ^ taken yet) is: 

511 = -{c/af - 2{c/a)u;oy + Y^Hn - 2{c/a)X^YHri + {X^yfH^^ 

9i2 = -^"12 + {X3y/Y)H23, 9i3 = DuqYHis + DX^yH^^, 

922 ^H22/Y\ g23^DH23/Y, g33^D^H33, Y^awoy + c. (3.4) 

Similarly as before, in the limit a;o the Hij will depend only on t, and the subse- 
quent limit of zero shear is Hu = —CiiR^{t) + (c/a)'^, other Hij{t) — —CijR'^{t), C33 = 1. 
Proceeding exactly as in sec. II, we then find that the hypersurfaces t = const will have 
constant curvature when A3 — > 0; the resulting limit is the Friedmann k = model, as ex- 
pected for type III. The hmits C3 — > and A3 — > should be tuned so that C3/ A3 > 0, 

wo— >0 

e.g. A3 = L^uq^/'^. 

The case 1.1.2.2 (eqs. (4.12) - (4.33) in Ref. 2) is of Bianchi type II. It is known from 
the paper by Ozsvath^. and from Theorem 3.1 by King and EUis^, that dust models of type 
II have zero rotation. However, that thesis was proven with use of the Einstein equations 
in Ref. 7 and of the Ellis evolution equations^ in Ref. 8, that inchidc consequences of the 
Einstein equations. In the approach of Refs. 1-3, the Einstein equations were not used. 
Moreover, the constant A3 plays the role of the tilt parameter here - with A3 = 0, the 
metric becomes stationary (the orbits of the symmetry group become timelike and tangent 
to the velocity field of the dust), and this case is not covered in Refs. 7 and 8. This is why 
the case 1.1.2.2 could show up in our consideration. This observation implies a warning: 
the existence of a Friedmann limit of the metric does not guarantee that the Einstein 
equations will allow a rotating generalization of a given Bianchi type of the corresponding 
Friedmann model. A rotating dust solution and the Friedmann solution may turn out to 
be two disjoint subclasses within that type. 

The metric (eq. (4.13) from Ref. 2) transformed back to the Plebaiiski coordinates (by 
the inverse of (4.4) from Ref. 2) is: 
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912 = hi2 + A3y/i23, S'13 = C3(/ii3 + Xsyhas) , 

{922: 923: ^ss) = (^22, C'a/ias, Cs^has), (3.5) 

where the /i^ are arbitrary functions of the T from (3.2). The hmit of zero rotation and 
zero tilt is achieved after the reparametrization: 

(y^Cs) = LOo(ij,D), 

{hl2: hi3) = {Hi2, Hi3)/uJ0: (/i22, /i23, ^Ss) = (-^22, H23, H33) / UJo^ . (3.6) 

and the reparametrized metric is: 

9n = i^oyf + hn + SAg^i^ig + {XsyfH^s 
912 = H12 + \zyH2z, 9i3 = D{Hi3 + XayH^s), 

{922: 923: 933) = (-f^22, M23, D'^Hsa). (3.7) 

The k — Priedmann hmit will result now when a;o = 0, Hij — —CijR^ and A3 = 0. 

The theorem by King and Ellis mentioned above implies that a;o = will follow when 
(3.7) is substituted in the Einstein equations. 

The Killing fields for the metric (3.1) are: 

= cxS^'o + axS'^i - {ay + c)S'^2 + {c\^/C^)x5''^. (3.8) 

(The Killing fields for (3.5) result when a = above.) After the reparametrization (3.3), 
in the limit a;o 0, the basis (3.8) becomes: 

^(1) = ^(3) (1/A3)^r3) ----- <^"3 

lm = -{^olc)k%) > ^2 - (A3/Z^)x5«3. (3.9) 

In the Friedmann limit A3 0, the generators (3.9) become a Bianchi type I algebra. 



IV. Cases 1.2 and 2 of Ref. 2. 

All of these allow both the /c = and the k = —1 Friedmann limits. 

Case 1.2.1.1 is of Bianchi type VI^ with the free parameter (6^ -|- P)/{h'^ — P) (there 
is a typo in Ref. 2). In this case (eqs. (5.6) - (5.7) in Ref. 2), the transformation back 
to the Plebanski coordinates is given by (5.5) from Ref. 2, with the roles of and x'"^ 
interchanged. The resulting metric is: 

9ii = hn - 2bt{y - bhu) - 2Zhi3 + {btf{b''h22 - 1) - 2bHZh23 + Z^hss: 

9l2 = hi2 + bHh22 - Zh23: 9l3 = hi3 + bHh23 - Zh33, 
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{922, 923, 933) = {h22, /i23, /iss), Z := fz, (4.1) 
where b and / are arbitrary constants, and hij are arbitrary functions of: 

T = t + y/b. (4.2) 

The hmit of zero rotation (that will automatically untilt the model) is achieved by the 
repar ametrization : 

y = ujoy, h22 = H22/i^o, h23 = H23/UJ0, 

hn = Hn + {bTf - 2b^Thi2 - b^T^H22/uJo^ 
hi2 = Hu/uJo - b^TH22/uJo'', his = i^ia - b^TH23/uJo, (4.3) 
and the repar ametrized metric is: 

511 = (oooyf + Hn - 2byHi2 - 2ZHi3 + {byfH22 + 2hyZH23 + ^'/laa, 

9i2 = H12 - by H 22 - ZH23, 913 = - by H 23 - Zh^s, 

{922, 923, 933) = {H22, H23, hss). (4.4) 

With cuq = 0, the shearfree limit will result when all hij = —CijR'^{t), and then the A; = — 1 
Friedmann model results when b = f ^ 0. The A: = Friedmann limit results when 
b = f = 0. This is the first instance where the coordinates of the k = —1 Friedmann limit 
come out rather exotic. From now on, this will be the rule. 
The KiUing fields for the metric (4.1) - (4.2) are: 

k(i) = S^u k^2) = e'"(5"o - W^), fcfs) = e^"5"3- (4-5) 

In the k = —1 Friedmann limit that will result by the first of (4.3) and b — f, k([-^ remains 
unchanged, /c^^ becomes e''^5"(3), while /c^^ is replaced by: 

= {-u^o/b)kf,^ __> e'^S^. (4.6) 

This is of Bianchi type V, and in the further limit 6 = / = it becomes type I. 

In the case 1.2.1.2 (eqs. (5.10) in Ref. 2), which is of type IV, the transformation back 
to the Plebahski coordinates is given by eq. (5.9) there. The whole further calculation is 
similar to (4.1) - (4.4) above. Instead of the last formula in (4.1) we have: 

Z ■.^ct + bz, (4.7) 

where c is one more arbitrary constant, and in (4.3) we have: 

hu = Hn + (bTf - 2b^Thu + 2cTHis - b^T^H22/uJo + {cTfh^^, 

hi2 = Hi2/u;o - b^TH22/u;o^ + cTZ/sa/c^o, = His - ^'TZ/sa/c^o + cT/iga- (4.8) 

The repar ametrized metric is: 

9n = {i^oyfi^ + c^hss/b'^) + 2cuoy{Hi3/b - ^i^aa - zh^s) 
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+Hu - 2byHi2 - 2bzHi3 + {byfH22 + 2b^yzH23 + {bzfhs, 

912 = ujQ{c/b)yH23 + Hi2 - byH22 - bzH23, 

gi3 = uJo{c/b)yh33 + H13 - byH23 - bzh^^, {922, 923, 933) = (^^22, -f^23, H33). (4.9) 

The limit a;o — > of (4.9) is the same as the hmit — > of (4.4) with b = f. Hence, 
the k — —1 Priedmann hmit will result from (4.9) when a;o = and hij — —CijR^{t), 
without any further limitations. The k — Q Priedmann limit will result when 6 = in 
addition. 

For the case 1.2.2.1 (eqs. (5.17) - (5.18) in Ref. 2), which is of Bianchi type VI/j, 
the subcase C = j + a = is identical to the subcase c = of case 1.2.2.2, and so the 
Priedmann limits will be the same (see below). 

The case 1.2.2.2 (eq. (5.19) in Ref. 2), which is of Bianchi type IV, allows the special 
case c = 0, where the Bianchi type becomes V. This special case was investigated in detail 
in Ref. 4, and it was shown there how the k = —1 Priedmann limit is obtained. In order 
to obtain the /c = Priedmann limit, one has to apply the following transformation and 
rescaling to eq. (3.5) in Ref. 4: 

y = e"", K = k/a, (4.10) 

and then take the limit a — > 0. 

All the subcases of case 2 in Ref. 2 have matter density constant along the dust flow: 
in case 2.1.1 (type I) and both cases 2.1.2 (types II and I), the velocity field is tangent 
to the symmetry orbits, in case 2.2 (type I), the velocity field coincides with one of the 
Killing fields. Therefore, no Priedmann limits will exist there. 

With this, all cases of Ref. 2 are exhausted. 

V. Case 1.1.1.1 of Ref. 3. 

In the cases considered in Ref. 3, each of the 3 Killing vectors is linearly independent of 
the velocity and rotation. However, the 5 vectors existing in each 4-dimensional tangent 
space to the manifold cannot form a linearly independent set. The 3-dimensional space 
spanned by the Killing vectors, K3, must intersect with the 2-dimensional plane spanned 
by the velocity and rotation, H2, along at least one direction. In the models of Ref. 1, the 
whole H2 plane was a subspace of the K3 space. In consequence, the velocity was a linear 
combination of the Killing vectors, and so those models were stationary. In the models of 
Ref. 2, considered up to now, the plane H2 and the space K3 intersected along the direction 
of the KiUing vector ^"3^ = + (As/n)^". Prom now on, the line of intersection will 
not coincide with any Killing direction. Hence, in each case an equation of the following 
form will have to hold: 

ai/c"i) + a2k'^2) + «3fc°3) = + ^2^", (5.1) 

where and bi are functions on the manifold. Note that if 62 = 0, then the velocity field 
is tangent to the symmetry orbits, and in consequence such a model has zero expansion 
and matter-density independent of the comoving time (the metric may depend on the time 
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only because in general the metric has shear). Hence, 62 is a measure of the tilt of the 
velocity field with respect to the orbits. Its relation to the tilt defined by King and EUis^ 
will be explained below (see after eq. (5.6)). 

The case 1.1.1.1, given by eqs. (2.28) - (2.29) in Ref. 3, is of Bianchi type VI/,. The 
transformation back to the Plebahski coordinates is given by eq. (2.27) in Ref. 3, and the 
result is: 



bp{b - f ) 



fl 



1 2 



bf3{b - f) 
fl 



33, 



23, 



912 = V/F + hi2 - Vh22 - ,^ ' ,M h 

bl3[b - f) 
fl 

913 = ^13 - Vh23 - 7-77 77^^33, 

b(3{b - f) 

922 = -l/b^ + h22, {923, 933) = (/^23, h3), (5.2) 
where b, /, (3 and 7 are arbitrary constants, the hij are arbitrary functions of the argument: 

T^t + y/b-P{b-f)z/r (5.3) 

and U and V stand for: 

U^bt + y, V^bft + {b + f)y. (5.4) 
The Killing fields for the metric (5.2) are: 



k%)=e'^{-d^ + bd^). (5.5) 



■^(3) 

From (1.8) it follows then that: 

be-f^k^2) + /e-'"^f3) ^ (b - f)u^ + b/{nP)]w^. (5.6) 

This is the equation (5.1) specified for the case 1.1.1.1. As remarked above, when 7 = 0, 

the velocity field becomes tangent to the symmetry orbits. (With 7 — > 0, the argument of 

hij given by (5.2) has to be redefined so that it becomes Z — 7T > — /3(6 — f)z.) 

7— *o 

This means that the parameter (7//3) is a measure of the tilt of the velocity field with 
respect to the symmetry orbits. Indeed, the measure of tilt defined by King and EUis^ is 
proportional to {l/P)- They defined the hyperbolic angle of tilt P by: 

cosh /3 = rt^riQ, (5.7) 

(the difference in sign from their paper is a consequence of the difference in signature), 
where is the unit vector normal to the orbits of symmetry. This definition of (3 makes 
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sense only when both and are timelike vectors; the cases of being null or spacelike 
are not considered in Rcf. 8. However, u'^Ua is a measure of the tilt also for nontimclike 
n". In particular, when is tangent to the symmetry orbits, u'^Ua = 0. The vector is 
related to the Kilhng fields by: 



where: 



In our case then: 



ria = Nj^-g^^N^^N- := Nj\\Nl (5.8) 

Na = --^£a/37<5^fl)%^(3)- (5-9) 

u^'ria = \\N\\-'N^u'' = __l_e(''+/)-7//3. (5.10) 

Analogs of (5.1) and (5.6) will exist in every case considered from now on. In the 
models of Ref. 2, considered up to now, where the Killing field k^^^ always had the form 

^(3) ~ C'a'^" + (^a/"^)""^"? -^3 was a measure of the tilt. 

For calculating the limit of zero rotation and zero shear, the following reparametrization 
is useful: 

(y, P) = uoiy, -B), /i22 = -f^22/^^o^ ^23 = -^^23/^^0, 
hu = {fTf + //n + 26/T/112 + 2{f/D)Ths/u;o 
-{hfTfH,,/u;^ - 2{hf/D)T^H,,/uj^ - {f / D^)T^h,,/uj,\ 
hi2 = H,2/u;o + bfTH^^/oJo'' + {f / D)TH2^/uJo\ 
hi3 = + bfTH^s/uo + {f / D)Thss/ouo, 

D:^f{h-f)Bh. (5.11) 

The reparametrized metric is: 

511 = -2ftWuJo + [{Dzf - 2{f/h)yWW + - 2hWH^2 
-2fzH,, + {bWfH22 + 2bfzWH23 + {fzfH,,, 

912 = Ut/h)uJo + [ih + f)/h'']yUo^ + i/i2 - hWH22 - fzH2s, 

gi3 = His - bWH23 - fzhss, 

{922: 923: 933) = {H22: H23, /t33 = -f^33), W := y + Dz. (5.12) 

The k — —1 Priedmann limit results now from (5.12) when uq = Q (after which all 
depend only on t), and: 

Hij^-djR'it), C33 = l, b^f. (5.13) 
(The last of (5.13) imphes D — 0.) The k = Priedmann hmit results when in addition: 

b = f = 0. (5.14) 

The reparametrization (5.11) transforms the Killing fields fc^-j and fc^^ from (5.5) as 
follows: 

Zf2) = (W7)fcr2) --> e^^[-(6//3/7)5"(2)+ 5^^(3)1 

/^3) = (c.o/6)A:^3) --> e^-5"(2). (5.15) 

In the k = —1 Friedmann limit {b = f ), together with k"-^-^, this becomes a Bianchi type V 
algebra, and in the k = Friedmann limit (6 = / = 0), it becomes type I. 
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VI. Case 1.1.1.2: of Ref. 3. 

This case is given by eqs. (3.7) - (3.11) in Ref. 3. It is of Bianchi type VIII or VIq (when 
g ^ OT g — 0, respectively), so only the k — Priedmann limit may exist here. The limit 
of zero rotation and zero tilt can be considered without transforming the metric back to 
the Plebahski form, but the 3 subcases have to be considered separately. 
The Killing fields in this case are: 

= e-^^''[2gyS^ + a,{gy' + 2B)5^ + 5^]. (6.1) 

The analog of eq. (5.1) is: 

k^2) - May^ + 25)e-2"i^A;^3) - 2gye-"''' kf^^ = e-°i^[4Sii" + {8c-fai/n)w% (6.2) 
and, consequently, the King-Ellis measure of tilt is: 

y/^N^u'^ = ai (6.3) 

{Na is given by (5.8)). 

The argument of the arbitrary functions in the metric is: 

T^t + y/a,--^z. (6.4) 

In case I {gB ^ 0), the reparametrization needed is: 

y = uoy, B = uq^^'^B, ai = LJo^^^ai, (6.5) 

(^13, ki3) = (Gi3, Ki3)uJo h22 = ^22/^^0^, (6.6) 

The full reparametrized metric, with cuq ^ 0, is rather complicated here, so only the limit 
a;o — > will be quoted: 

ds^ = dt^ + /iiidx^ + G22^y^ - -^dydz + h33dz\ (6.7) 

2Bai 

The A; = Friedmann limit results when further hu = G22 = h^s = —R^{t). 

As seen from (6.1), the symmetry group becomes Bianchi type I in the limit a;o — > 
after the reparametrization (6.4) (the Killing field /c(3) has to be replaced by Z(3) = ujQk(^3) 
in order that the limit is nonsingular) . 

In case II {B — 0), eq. (6.5) remains unchanged, while (6.6) is replaced by: 

ihi2, h23) = (G12, G'23)a;o /i22 = G22/<^o'- (6.8) 
The limit a;o ^ of the reparametrized metric is: 

ds^ = dt^ + hiidx^ + 2Gi2dxdy + 2hi3dxdz + G'22dy^ + h^^dz^ , (6.9) 
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where all the metric components depend only on t. The A; = Priedmann limit is here 

Gi2 = hi3 = 0, hn = G22 = = -R\t). (6.10) 

In case III (g = 0, Bianchi type VIq), the k — Friedmann limit results again by (6.5), 
(6.8), (6.9) and (6.10). 

Note two typos in Sec. 3 of Ref. 3: in (3.7) the correct formula for is: 

and in (3.10), the correct formula for g23 is: 

923 = -2gzhi3 + h23 - axgz^h33. (6-12) 

VII. Case 1.1.2.1 of Ref. 3. 

This model is of Bianchi type VII/j, and is given by eqs. (4.19) - (4.23) in Ref. 3. Two 
formulae in (4.23) had typos, the correct expressions are: 

= e^^+M''[Wh^2 - {l/D){h + /) cos(L'x/2)/ii3], 

= c(^+^)^{[7'(6 + ffh33/D^ + 1] COfi\Dx/2) 

-2{j/D) {b + /) cos{Dx/2)Wh23 + W^h22] (7.1) 

The transformation back to the Plebahski coordinates is given by (4.21) in Ref. 3, and the 
resulting metric is: 

911 = + /ill + 2C//112 - 2(7/^) (6 + /)y/ii3 + U''h22 

-2(7/D)(6 + f)yUh23 + {{llD){h + /)y]'/i33, 
912 = -2/112 + 2(7//^)/ii3 - 2?7/i22 + 2(n/D)Uh23 
+2{j/D) {b + f)yh23 - 2{j/Df{b + f)yh33, 

9l3 = hi3 + Uh23 - {l/D){b + /)y/l33, 
g22 = 4/l22 - 8(7/L>)/l23 + 4(7/L>)2/l33 , 923 = -2/l23 + 2(7/L>)/l33, 

to = /i33, ^ := \[{h + /)' + + 2(6 + /)y, (7.2) 
where b, f D and 7 are arbitrary constants, and hij are arbitrary functions of the variable: 

The Killing fields for the metric (7.2) are: 

kfi) = A;f2) = e(^+^)^/2[cos(Da;/2)5o" - ^W-52" - 7sm(I^a;/2)53"], 
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W ■.= {h + f) cos{Dx/2) - D sm{Dx/2), 

V cos{Dx/2) + {b + f) sm{Dx/2). (7.4) 
The analog of (5.1) is here: 

Vkf^) - Wk'(^^ = r>e(^+^)^/'ii° - ^(6 + /)e(^+^)^/'«;'*, (7.5) 

and the King-EUis measure of the tilt is: 

V^W^N^ --lib + f)W'^f^^- (7.6) 

Eq. (7.5) shows that with 7(6 + /) = 0, the model should be nonexpanding. This is so 
indeed, but in order to be able to consider the subcase 7(6 + /) — > 0, we have to take 
7(& + f)T as the argument of hij in (7.2) instead of the T given by (7.3). 
We define: 

E := (b + fy + D\ (7.7) 
and then the reparametrization needed for the limit of zero rotation and zero tilt is: 

(y, D) = ujo{y, d), /i22 = i?22/t^o^, ^23 = -^^23/^0, 
hn = Hu - ETH^^/oJo + [l^T^ ^22^', 

hi2 = H12/UJQ — -ETH22/i^Q^ 1 hi3 = Hi3 — -ETH23/UJ0. (7.8) 
After the reparametrization we have: 

Si,2-=b + f + 6i^2{uJody/{b + f) > 6 + /, £i = +l,£2 = -l, 

OIQ— >0 

E := S2y - ^{d/^)Siz, 

911 = {c^oyf + Hn + 2J:H,2 - 2{^/d){b + f)yH^^ + Y?H22 
-2(7/d)(6 + /)yE//23 + [{l/d){b + f)yfh^z, 
gi2 = -2i/i2 + 2{-f/d)Hi3 - 2E//22 + 2{-f/d)yH23 

+2(7/rf)(6 + / + S)i/23 - 2(7M'(& + f)yhs, 
gi3 = i/i3 + yT.H23 - il/d){b + f)yh33, 

922 = ^H22 - Kl/d)H23 + 4(7/(i)2/i33, 
5^23 = -2-ff23 + 2(7/ci)/i33, 5^33 = h^z := i?33- (7.9) 

In the limit tuo — > 0, all the Hij will depend only on t. The limit of zero shear is then 
obtained by: 

H,j = -Ci,R^it), C33 = l. (7.10) 
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To obtain the Friedmann limits, a further reparametrization of the constants Cij is neces- 
sary. We define: 

C23 = -D23 + 7/ d, D22 = C22 — C23 , 
D12 = (C12 — C13C2Z) I D22-I Dii' = Cii — Ci3^ — Di-^. (7-11) 
The metric (7.9) may then be written: 

^ 2 



ds^ = dt^ - {DnRdxf - 



-Di2-{b + f)D22 [y-^^ 



dx + 2D22dy 



Ci3 + {b + f)D2sy -^{b + f)(l + ^023) z 



dx - 2D23dy + dz 



(7.12) 



The k — —1 Friedmann hmit results now when d — 0, the k — limit results when b+f — 
in addition. The first of (7.11) was necessary to eliminate 7/d from (7.9) so that the hmit 
d could be subsequently taken. 

We have found above (after eq. (7.6)) that 6 + / = corresponds to zero expansion. 
This is so when 6 + / — > with other parameters uchanged. In considering the k = 
Friedmann limit, 6 + / is set to zero after the limit d — > had already been taken. In order 
to make these two limits compatible, we have to assume that 6 + / — > slowly enough so 
that D/{b + f) ^ and y/{b + f) ^0. With the reparametrization (7.8), this is achieved 
when b + f = Buoif ^ where < 5 < 1. 

After the reparametrization (7.8), in the limit ujq — >■ 0, the Killing fields become: 



(1)' 



l'(2)^-[2u;o/{b + f)]kf2^ > e(''+^)-/2(5"(2) 

a;o— >u 



^^3) = (l/7)^r3) 



-> e 



{b+f)x/2 



d 



l + -{b + f)x 



(5" (2) + (5" 



(3) 



(7.13) 



This is a Bianchi type IV algebra, and in the Friedmann limits k — —1 {d = 0) and k — 
(d — b + f — O) it becomes type V and I, respectively. 



VIII. Case 1.1.2.2 of Ref. 3, Bianchi type IX subcase. 

The case 1.1.2.2 contains three different subcases that are of Bianchi types IX, VIII and 
VIIq. The type IX subcase requires some adaptation of the formulae given in Ref. 3. 

For type IX, g/c > 0. Then, as seen from eq. (5.16) in Ref. 3, B/c < 0, or else (5.16) 
would lead to a contradiction. These two inequhties imply that gB < 0, while eqs. (5.26) 
and (5.27) in Ref. 3 are adapted to the case gB > 0. Hence, a re-adaptation of these 
formulae to type IX is necessary first. We define: 

B := -B, A = iX, ki2 = iki2, k23 = ik23 (8.1) 

(the overbars simply denote new symbols that will be real), so that instead of (5.16), (5.23), 
(5.26) and (5.27) from Ref. 3, we obtain: 
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R^2cD%/{BK), 



K'^Rdy 



v^Bt + 2cD-fz, U = hi2 sin(2A't;) + ^008(21^); 



(8.2) 



911 = + K^Hn + A^yKH 



B 



13 



gB 



H- 



33) 



D 

9i2 = Hi2 + 2=—yH23, 



^13 = KH,3 + 2^|/i733, 



H 



BK'' ' " B 
g22 — H22/K'^, g22, — H2z/K, Qzz — H^z] 



(8.3) 



11 



25^X 

Hi3 



U + /ill, -f^i2 = /ii2 cos(2Ai;) — /ci2 sin(2Ai;) 



252A 



[h23 sin(At') + k23 cos(A^;)], 

2 



22 



U + 



gB iicY 



:hn + 



CD^ 2c2D6— ■ 5/^2 

-f^23 = h23 COs(Xv) - k23 Sm{\v) , H33 



h33, 



h 



33; 



U) 



where B, c, D, g and 7 are arbitrary constants, and all the hij, ku and /C23 are arbitrary 
functions of the argument: 



B 



2cD-f 



z. 



(8.5) 



Eqs. (8.2) - (8.5) are written in the Plebanski coordinates. 
The KiUing fields for the metric (8.3) - (8.4) are: 



i^il) — Ol 1 



7 



k^^ = cos(Dx/2) {K - yK,y + K,, 6^ + 



A;^3) = sin(Dx/2) 



{K-yK,y)5^ + K,y5? + 



7 



DK 



+ -DKsm{Dx/2)S^, 

1 



- -DKcos{Dx/2)5^. 



(8.6) 



(Note: the first commutator in eq. (5.20) in Ref. 3 should have a minus on the right-hand 
side.) The analog of (5.1) here is: 

cos(L>a;/2)A;^2) + sin(L>a;/2)A;^3) - K,y A;^^ ^ (K - yK,y )u'^ + +[^/{DKn)]w". (8.7) 

In agreement with this, the King-EUis measure of tilt is here: 

V^TV^ii" = 7/2. (8.8) 

Note that: 



K 



gy 



K-yK„ 



B 



(8.9) 



icD^K' " ^"'^ 2cD^K' 
The case presently considered is the only one of type IX in the whole classification. 
Therefore: 
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1. This is the only place where the k = +1 Friedmann model will appear as a limit. 

2. The models represented by eqs. (8.2) - (8.6) include those considered by Godel^'^. 
(Ours are in fact more general because the tilt of the symmetry orbits with respect to the 
velocity field is an arbitrary parameter here.) We shall deal with this point further on. 

For later considerations, it will be convenient to reparametrize the metric (8.4) once 
more, as follows: 

G := (/ii2' + k^y^^ , /ii2 = -Gsin(2/3), ^ = Gcos(2/3), 

F := (h23^ + k^''y^\ h23 = Fcos(a), ^=Fsin(a), (8.10) 

where G, F, a and P are new functions of the T given by (8.3); and also to transform the 
coordinate y by: 



y = pB/gcos{^). (8.11) 

From now on, the x^-coordinate will be i). Then, in order to set the rotation and the tilt 
to zero, the following further reparametrization is needed: 

B^W, 7 = ^^o, G^G/c^o, /ill = Gn/c^o'- (8.12) 
Let us note that: 

A^; = — ^ > — ^ — z. (8.13 

2x/2[(cuo&)2 + (2cD/i)2] 4V2c/i ^ ^ 

The metric (8.5), reparametrized by (8.10) - (8.12), becomes: 

2^ 2 2 26 ^ . 2 ,chD^ /2c ^ . ,t n • o 

011 = — COS vcuo + , ^o Giism i7 — 4 W — FsmfAi' + a) smvcosv 

g AcD^ g \ b 

-dJ—Gcos(2Xv + 2(3) sin^ d + S^^^^^hss, 
V 2«/ bg 

gi2 = ^j2b/gGs\n{2Xv + 2/3) sini? - 4:V2{chD^ / g)^/^F cos(Xv + a) cosi? 

chD [c — 

gi3 — 2\/2—=hs3 COS'!? — D^J-Fsm(Xv + a) siwd, 

Vbg V 9 



922 



Gil + 32^^/i33 + 2^2^ J -a cos(2A^; + 2/3), 



cD^ bg D\ g 

g23 = -2D^Fcos(Xv + a), ^33 = hs^. (8.14) 
In the limit uuo ^ 0, the velocity field = will have zero shear when: 

a — const, (3 — const, 
(G,F,G'ii,/i33) = -(Gi2,G23,Gii,l)i?2(t). (8.15) 
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With use of (8.13) it may be verified now that the k — +1 Priedmann hmit will result from 
(8.14) - (8.15) when ujq ^ and: 

Ci2 = C23 = Cn=0. (8.16) 

The resulting representation of the Friedmann model (again an exotic one) is identical, up 
to rescalings of coordinates, to the one derived by Behr^^ 



ds" = dt^ - R\t) 



8 ^ . ^ dx^ + 4V2^= cosi?dxd^ + 32^-^d'&^ + dz^ 



(8.17) 



bg y/bg bg 

The A; = Priedmann model will result from this after the transformation/reparametrization: 

•d = aiccos{ky), h = H/k, x = kx', (8.18) 

in the limit A; ^ 0. 

For the k = +1 Friedmann limit, the algebra of the Killing fields (8.6), suitably trans- 
formed by use of (8.9), (8.11) and (8.12), is still of Bianchi type IX. For the k — limit, 
the algebra {/(i), /(2), ^(3)} := A;{/c(i), /c(2), /c(3)} is of Bianchi type 1 when A; — >• 0. 

As stated above, the class of models defined by (8.2) - (8.5) must contain the one 
considered by Godel in Ref. 10. This is so because two of Godel's assumptions (dust 
source and nonzero rotation) place his class within our collection, and the third assumption 
(compact spaces t = const, i.e. Bianchi type IX; the Bianchi classification and terminology 
had not yet been in common use in Godel's time) uniquely points to the subcase I of our 
case 1.1.2.2. Godel presented several properties of these models in the form of theorems, 
but mostly without proofs and almost without formulae. It would be an interesting exercise 
to see how Godel's theorems apply to the explicitly given metric (8.2) - (8.5). 

In particular, one of his statements seems to need a refinement. He said that there 
exist 00^ rotating solutions satisfying all his requirements. This means that the collection 
of all solutions of the Einstein equations for (8.2) - (8.5) should be labeled by 8 arbitrary 
constants. One can understand how this happens from Ref. 4, where the Einstein equations 
were investigated for an equally general Bianchi type V class. Of the 6 unknown functions 
in the initial metric, one (^^33 in Ref. 4) is determined by an algebraic relation, two of 
the Einstein equations are of first order and can be used to eliminate two more functions, 
and then the remaining 3 functions obey equations of second order. This gives 8 constants 
indeed. However, the tilt parameter (7 in (8.7) - (8.8)) is one more arbitrary constant that 
is contained in the metric even before the Einstein equations are considered. 

Rotating dust models of Bianchi type IX were considered by Behr with simplifying 
assumptions about the metric. Similarly as in Ref. 4, the main conclusion seems to be 
that whatever one does with the Einstein equations, no solution comes within sight. 

IX. Case 1.1.2.2 of Ref. 3, Bianchi types VIII and VIIq. 

The subcase of case 1.1.2.2 that corresponds to the Bianchi type VIII is defined by: 

g/c<0 (9.1) 
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in eqs. (5.16) - (5.27) in Rcf. 3. Then, B/c and, consequently, Bg can have any sign at 
this point. Only the A; = Friedmann model can be contained as a subcase here. 

The cases Bg ^ and Bg = have to be considered separately. When Bg ^ 0, we 
take eqs. (5.23) - (5.28) in Ref. 3 with the following specializations: 

hi2 — ki2 — /i23 — ^23 — 0, 

B^bujo^, {'j,y) = {h,y)LJo, /in = GiiM^- (9.2) 

Then: 



~ ~ 1 / gr + 2b 
K = uoK, K=—^ 

bg 

912 = 0, gi3 = -2{cDh/b)yh33, 

922 = ( ~ 2^2/^6 ^11 ~ ^bD^^^^ J ' ^ °' " ^^'^^ 

It is now seen that the proper signature will result only when: 

b/c > 0. (9.4) 

In order to obtain the k — Friedmann model from (9.3), we then rescale the constants 
again as follows: 

b^boVg, D^dg^'\ h^Hg^l\ (9.5) 
and take the limit — > 0. The limiting metric is: 

d.^ = de + Gn (2^d^' + J,^f) + h3,dz' (9.6) 

and in the hmit of zero shear, Gn = Cn/isa = —CiiB?{t), this becomes the A; = 
Friedmann models indeed. 

The Killing fields for this case are given by (8.6). After the rescalings (9.2) and (9.5), 
the following basis of the symmetry algebra is obtained: 



This is of Bianchi type VIIq. 

When Bg — 0, and the Bianchi type is VIII, we must have: 

B = 0^g. (9.8) 
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The metric is then found from eqs. (5.28) and (5.23) - (5.26) in Ref. 3, suitably adapted. 
With B = 0, the arbitrary functions depend only on t. The metric needs then to be 
rescaled as follows: 

(7, y) = {h, y)uJo, (/ii2, hs, /laa) = {G12, G13, G2^)/ujq, hu = Gh/uq (9.9) 
and the result of the rescaling is: 

gii = i^oyf - ^^'^11 + ^^^^^12 - 2^~c/ghGi3 



^cDh , Dg o, 
9is - i^^-9lcyG,, + ^y zG,, + + ^yz 
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^22 = -4^/122 + 2{z/y)G23 + {z/yfh. 
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5-23 = G23 + (^/y)/i33, 933 = h3- (9.10) 

In the limit c<Jo ^ one term in gu disappears and the hij depend only on t. The 
shearfree limit is then attained when: 

Gij^-GijR\t), h33^-R\t), h22^-G22R\t). (9.11) 

To find the Priedmann limit we then assume that: 

Ci2 = Ci3 = C23 = 0, h^HD, g^-cG^D\ (9.12) 

where H and G are new constants (the last definition takes into acccount that g^/c < in 
type VIII), and let a;o — 0, -D — > 0. The resulting metric is: 

ds''^de-Gii{^GyR^\x'-G22{2^ dy' - i?^ |^-2^dx + |dy + d^^ .(9.13) 

The A; = Priedmann limit results from this when: 

y = e^ C22 = {D22/k)\ k^O, (9.14) 
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With (9.8), (9.9) and (9.12), the KiUing fields become: 

— 0^ , 

A;" = -Gcos{Dx/2)5'^ + \GDysin{Dx/2)5^ + ^ cos(L'x/2)5", 
^ ' 2 4 Gy 

1 1 1H 

k^^^ = -Gsm{Dx/2)6'^ - -GDy cos{Dx/2)S^ + —sm{Dx/2)S^. (9.15) 

Before the hmit D — > can be taken, k(^s) needs to be redefined by /c^g) = {1/D)k(^y The 
basis in the hmit becomes: 



k'(,) = 5i. kt2) = 7;GS'^ + -^Sl lh=^2 (9.18) 



1 1 Hr 

k[% = -Gx5? - -Gm + -^51 (9.16) 

This is of type VIq. We now transform y by (9.10), and redefine /c'^g^ once more: 

11^ = -{Ak/G)k'^,y (9.17) 
In the hmit A; — > 0, the foUowing basis then results: 

{u being now the x^), which is clearly of Bianchi type I. 

Finally, when the Bianchi type is VIIq (i.e. = 0), the metric results by a simple 
specialization of eqs. (5.23) - (5.25) and (5.28) in Ref. 3. In this case necessarily b/c < 
and: 

= —W— = const. (9.19) 
D V 2c ^ ' 

The rescaling that will allow to calculate the limit a;o ^ is: 

y) = (&, y)<^o^ D = duo, 

{hi2, /123) = (^12, 6-23)^, h22 = G'22/u;o'. (9.20) 
The argument of the arbitrary functions must then be redefined so that it becomes: 

u = u/(2cD'y) = t + r——z > t. (9.21) 

2cD^ wo->o 

The limit a; — > of the metric is then: 

h 2 l—b cdP 

ds^ = dt^ - — — /iiidx^ - 2G'i2da;dy + -W —hi^dxdz - 2—C22df + 2G23dydz + hsadz^. 
2cd^ a V 2c b 

(9.22) 

The k = Friedmann limit results from here when shear is set to zero, i.e. when g^j = 
—CijR^{t). 

The basis of the KiUing fields in the limit 6i;o ^ is found as follows: 

%=*r, ir„=Hm (f3)=lim(-Jf-o*3")=*?. (9.23) 
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X. Cases 2.1 of Ref. 3. 

In the case 2.1.1 the transformation back to the Plebahski coordinates is the inverse of 
(7.16) in Ref. 3, and when applied to (7.18) there, it gives the following metric: 

911 = y'-{y + vf + /ill + (& + my + y)h^2 + {ilc)Vh^^\ 

+ |(^ + f)\{y + yfh^^ + 2(7/c)^(?/ + V)/l23 + (7V^/c)'/i33], 

512 = ^(2/ + V^) - /il2 - \{h + my + V)/l22 + (7/c)V/i23], 

513 = /il3 + ^(fe + my + V^)/i23 + (7/c)V^/i33], 
522 = -4/(6 + jf + /i22, 523 = -/i23 , 533 = /i33, (10-1) 

y :=^(6 + /)i + y, (10.2) 



where 



and the arbitrary functions depend on 

(Note two typos in Ref. 3: in eq. (7.18), the coefficient of Wh23 in (722 is 27, not 2b, and in 
(7.17), in the formula for m", there should be a (W/c) in front of S°'q.) The Killing fields 
for this metric are: 

k^i) = S^, k^2) = |cx(5o" _ c[l + 1(6 + /)x]52" + 7<^3"} 

kf,^^e('+f>/'[5^-^-{b + f)5^], (10.4) 

and they form a Bianchi type IV algebra. 
The analog of (5.5) is 

% - ^[1 + ^(6 + f)x]kfs) = e^^+ZW^ (^^«" + , (10.5) 

and the King-Ellis measure of tilt is 

x/=5«"iVa = -lib + /)7e('+^)^ (10.6) 
The redefinitions needed to make the limit a;o ^ finite are 

{y,c) =(jJo{y,C) 

hn ^Hn-^{b + ffT[hu + h/c)hs] 
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(b + ffT'[H22luj^' + 2(7/c)//23M + (7/c)'/i33], 



hi2 = i/i2/a;o - -^{h + ffT[H22lujo' + (7/c)//23/a;o], 

/ll3 = i^l3 - + ffT[H2^luJo + (7/c)/i33], 

h22 — H22/C0o^, h23 — H23/u!o, 



and the resulting metric is 



-(b + f)t + 2u;oy 
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(10.7) 



+{b + f)[YHi2 - zHn] + + ff[Y''H22 - 2YzH23 + z'^Hss], 
912 = ouot + - H,2 - lib + f)[YH22 - zH2s], 

913 = H^s + lib + f)[YH23 - zH^s], 

922 = H22 - 4a;o^/ ib + fY, 923 = -H23, 933 = ^33 = -f^33, 

Y:=y- Czh. 

In the hmit 0, aU H^^ become functions of t, and the Kilhng fields become 

— , 

1%, = lim k%, = e(^+^)^/' |-C[1 + + + T^^?} 

2 



(10.8) 



^?3^ = hm ^ ^ 

still of type IV. 

The shearfree limit of (10.8) is 



^oA:f3)) = e^'^M^^<^2. 



(10.9) 



a^i = -CnR'^it) + 



-I 2 



:ib + f)t 



other //i^- = -CijR^it), C33 = 1. (10.10) 

The k — —1 Priedmann model will then result when C = (and, consequently, Y — y), 
the /c = Friedmann model will result when b + f = 0, with no condition on C. Both 
limits can be easily taken also in the Killing fields (10.9), with C = they become of type 
V, with b + f = they become of type I. 

The case 2.1.2 was shown in Ref. 3 to be included in 2.1.1 as a subcase. 
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XI. Case 2.2.1.1 of Ref. 3. 

This case includes two subcases, Ay^O and ^ = 0, given by eqs. (9.11) - (9.15) in Ref. 3. 
Both are of Bianchi type VIII. The coordinates used there are those of Plebahski. 

With A 0, cqs. (9.14) in Ref. 3 are adapted to the case ^ < 0. However, when 
^ < 0, the hmit of constant curvature in the spaces t = const has a wrong signature. 
Therefore, the formulae must be re-adapted to ^ > 0. This is the result: 

gi2 = hi2 cos(2At') - ki2 sin(2At') + — [/i23 cos(Av) - sin(At')] , 

923 = [^23 cos(A^;) - k23 sin(Av)] , 5^33 = if33 = h^s, 
U :— /ii2sin(2Av) + A;i2 cos(2Av), V :— h23sm{Xv) + A;23 cos(Av), 

A^ = 2^/(4^2 + 7^)2, v = 2At + -fz, (11.1) 
and the hij are arbitrary functions of 

T^t- 2Az/-f. (11.2) 

The Killing fields for (11.1) are 

k^i) = s^, k^2) = (2^/y)<^o + i-^/y' + - ^ys^ + (i/yK, 

k'(,)=x6^-y5^, (11.3) 

The analog of (5.1) is: 

(^/y2 + xV2)A;[l) + k^^^ xk^^^ = {2Aly)u'^ + ^w'^, (11.4) 

yn 

and the King-EUis measure of tilt is: 

V^M"iVa = 7- (11-5) 

The redefinitions needed to make the limit ujq — > of (11.1) finite are: 

X 

{y,-f)^uJo{y,h), A=-{auJof, 

hn = Hn/uJo\ (/ii2, ki2) = (i^i2, Ku)/uJo. (11-6) 
Note that with (11.6) we have 

Xv > az/h. (ll-"^) 

Wo— >0 
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The reparametrized metric is 

gu = f [2hV/a^ + + a'^ [i/i2 sin(2Aw) + K12 cos(2Aw)]} , 

gi2 = {h/a^) [/i23 cos(A'u) - A;23 sin(A'u)] + H12 cos{2\v) - Ku sin(2A'u), 

913 = {y/a){V + hhss/a), 
922 = {-a [i/12 sin(2At;) + K^^ cos(2At;)] + a^iJn - (V«)'^33} , 

5'23 = [^23 COs(At') - A;23 sin(At')] , g2,z = -ff33 = ^33- 



(11.8) 



In the hmit o^o 0, all the arbitrary functions will depend only on t. 

The A; = Friedmann limit follows from (11.8) when the following further specialization 
and transformation is made: 



/i33 = -R\t), Hn = -(Cn + h^)R\t)/a\ 



H12 

and then the limit a - 



: K12 = h23 = k23 = 0, 

is taken. The metric becomes then 



Cn {dx'^ + du"^) + {gdx + dzf 



;il.9) 



(11.10) 



which is clearly the k = Friedmann model. 

The rescaling (11.6), followed by a;o — > 0, and the rescaling (11.9), followed by a — > 
transform the Killing fields (11.3) into an almost-standard Bianchi type I basis {k^^^ has 
to be replaced by l^^^ — a^k^^^, and k^^^ has to be replaced by = a/c^^^ before taking the 
limit a ^ 0). 

The case ^ = is given by eqs. (9-11) and (9.15) in Ref. 3. The rescalings needed 
there are 

y = ujoy, hii = Gn/ujo^, {hn, hia) = {Gu, Gxz)/uq. (H-H) 

The arbitrary functions hij depend only on t from the beginning. The limit a;o — > of 
the rescaled metric is 

ds^ = df + y'^Gndx^ + 2(^Gi3 + Gi2)dxdy + 2yGi3dxdz 

+y~^(^22 + 2zh22, + ^^^33)dy^ + 2y~^(/i23 + zh33)dydz + h^^dz"^. 
The A; = Friedmann model results now when 

G12 — Gi3 — 0, Gil — C'ii-R^(i), 



(11.12) 



hij — —CijE,"^, C33 — 1, C22 — 1/0^, 



0. 



(11.13) 

The Killing fields need not be reconsidered because ^ = is an allowed subcase for (11.3). 
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XII. Cases 2.2.1.2 of Ref. 3. 

In considering these cases, we first have to correct two errors. The first error is that the 
arbitrary constant yo actually must be equal to zero in all the formulae. The second error 
is that one subcase was overlooked - it needs special treatment and is not included in the 
formulae given in sec. X of Ref. 3. This special case is defined by 

^ = (12.1) 

and consequently /ii = and = j- It is because of /Xi = that some of the formulae do 
not apply to this case. 
The conclusion that 

a = 0, c=l (12.2) 

can be achieved by a change of the basis of the Killing fields is still valid. With (12.1) and 
(12.2), the solutions of eqs. (10.2) and (10.3) in Ref. 3 are 

P = -jy + M, L3 = 7, (12.3) 

where j, M and 7 are arbitrary constants. The resulting Killing fields are (by (10.6) from 
Ref. 3): 

^fi) = S?, fcf2) = MxS^ - 3x5^ + {jy - M)5^ + ^xS^ 

kf,) = M6S-j6'^ + ^6l (12.4) 

and they form a Bianchi type III algebra. 

The coordinates are still those of Plebanski at this point, so u°' and w°' and goa have 
their standard forms. The solution of the Killing equations is 

gu = m/j)y - {M/j f + ijy - Mfh, + 2{^/j)ijy - M)his + (7/jf /^as, 

gi2 = hi2 + -i^^^^^y 913 = {jy - M)hi3 + (7/j>33, 

922 ^ h22/{jy - Mf, g23^ h23/{jy - M), 5-33 = ^33, (12.5) 
where the hij are arbitrary functions of the argument 

T^t- {M/-f)z. (12.6) 

The Killing fieds (12.4) are a subcase of the general expression that will apply to the 
whole case 2.2.1.2 collection. The analog of (5.1) will be given further on for the whole 
class. 

The rescalings needed to find the nonrotating limit of (12.5) - (12.6) are 

iV: M) = uJo{y, m), hn = Hh/ujq^, 

(hu,his)^(Hu,H,s)/u;o. (12.7) 

The rescaled metric is 

gu = {2m/i)uJo^y - {mujo/jf + {jy - mfHu + 2{^/j){jy - m)Hi^ + {-i/jfh^^, 
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gi2 = Hi2 + -1— r, S'la = {jy - m)Hi2, + (7/^)^33, 

][]y-m) 

922 ^ H22/{jy -mf, 923 ^ h23/{jy -m), 933^ h^z, (12.8) 

The A; = Fricdmann limit is now obtained from (12.8) when shear is set to zero 
{hij — —CijB^it), C33 = 1), and in addition 

7 = hj, j ^ 0. (12.9) 

In order to make the hmits ujq ^ and j — > compatible, it has to be assumed that 
j oc uq", where < a < 1, e.g. a — 1/2. The Killing fields (12.4) become then an 
almost-standard Bianchi type I basis in the limit ujq — > 0, but k^^) to be replaced by 

70 _ , , -1/2 ua 
t(3) - 1^0 ft(3)- 

The case 2.2.1.2 consists of 3 subcases, each of a different Bianchi type. The subcase 
9 < is of Bianchi type VI/^, with the free parameter j/(j^ — AgY^'^. However, the 

paramctrization of the metric used in Rcf. 3 is inconvenient for calculating the Friedmann 
limit. It will be more convenient to rewrite it in the parametrization in which //i and //2 
appear symmetrically. Therefore, instead of (10.15) - (10.17) from Ref. 3, we will use the 
following formulae: 

U := Mcosh(L>y) + A^sinh(L>y), V := Msinh(L>y) + A^cosh(L>y), 
D (7-2/4 p.^ 7 

P = e-^y/^U, y = -^P - -e-^^/V, 

2^ 9 

9n = + hiiP^ + 2re-^^UVhn + {^V / D)e-^'' h^^, 
912 = hu + nV/U)h23, 913 = h^P + re-^'^/V/i33, 

922 = h22/P^, 923 = h23/P, 6-33 = ^^33, (12.10) 

where j, 9, M, N and 7 are arbitrary constants, y is one of the coordinates, Y is just 
a parameter used to represent P{y), and hij are arbitrary functions of the coordinate z. 
Since the coordinates used in (12.10) are not those of Plebahski, the u°' and 90a do not 
have their standard forms, they are given by eqs. (10.9) in Ref. 3. The transformation 
back to the Plebahski coordinates is the inverse of (10.8) in Ref. 3, after which we obtain 

912 = 9thu + /ii2 + rgt{V/U - \j/D)h,s + T{V/U)h23 - T^^t (l + ^) /133, 

922 = P-^ {{9tfhu + 2gth^2 - {jTg/D)t{gthi3 + /^23) + /i22 + giTgt/Dfhs} , 

923 = (l/P) [9thi3 + h23 - ^^33} , (12.11) 
with gii, gfi3 and 5^33 being the same as in (12.10). The hij depend now on 

T^t+^z. (12.12) 
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The Killing fields corresponding to (12.11) - (12.12) (and to all the other subcases of case 
(2.2.1.2)) are 



A;'(l) = 5?, = x(P - yP,y + xP,y 5^ - P5^ + xi^ / P)e-^'' 5^ 

= {P- yP,y + P,y 5^ + (7/P)e-^-^53" 



The analog of (5.1) (again valid for all the subcases) is 



kl) - P,y = {P- yP,y + 
and the King-EUis measure of tilt is: 

^^u^Na = -76-^'^. 
The rescalings needed to make the limit a;o ^ finite are 

(M, N, 7) = (m, n, h)ujQ, hu = Hu/uo^, hi3 = Hi^/ujq. 
In consequence of this we have 

{y,P,U,V)^{y,P,U,V)uo, 



;i2.13) 



;i2.i4) 



;i2.i5) 



(12.16) 



(12.17) 



where the symbols with a tilde are obtained from those on the left by replacing {m,N) — > 
{■m,n), and they do not depend on cuq. Also, from now on Y will be used as the x^- 
coordinate in place of y, so 

dy^cuoPdY. (12.18) 

The rescalings (12.16) have to be accompanied by the following redefinitions of other 
functions in the metric: 



h22 = -{gTfHn/uJo' - 2gThu + ^^22 + ^ {gT^H^^/uo + T/123) - giTgT/Dfh^, 

(12.19) 



^23 = -gTHis/ujo + H23 + ^^^^^33- 



The metric resulting after the redefinitions and the coordinate transformation is 



gu = M'' + H^iP^ + fe-^^^ [2UVH^3 - (h/D)hss 
gi2 = P {h,2 - {D/r)zHn + {j/2 - DV/U)zH 

5i3 = i^i3^ + fe-^''/2\//i33, 



13 
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g-22 = {Dz/rfHn - {D/t){2zH^2 + jz'H^z) + H22 + jzH^z + gz^h^^, 

923 = -{D/t)zH^^ + //23 + U/'2)zhss, ^33 = /iaa = ^33, (12.20) 

In the limit cuq — > 0, all the Hij will depend only on t. 

The limit of zero shear is then, as usual, Hij = —CijB?{t), C33 = 1, and the k — —1 
Priedmann limit results from (12.20) when, in addition 

h^HD, Ci3 = 0, D^O. (12.21) 

With D = we have U = m, V — n. The, again rather exotic, representation of the 
limiting Friedmann model is 

ds^ = dt^ - R^(t) (mDne-^^/^dx + DudYf - (D22RdYf 



Hn 

^-jY/2^^.^ + (C23 + j^/2)dy + dz 



w? — -n? 



(12.22) 



where 



m? — -n? 



1)22^ '■ — C22 ~ (^23^ ~ Di'^ . (12.23) 

The k — Q Priedmann limit results from (12.22) when j — 0. 

The rescaling (12.16) and the limit a;o ^ transform the KiUing fields as follows 

kl) = k^2) = (-j72 + DV/ij)x5'^ - 6^ + {h/U)xe-^''/^S^, 

kfs) = (-j72 + DV/U)5^ + {h/U)e-^'"^5t (12.24) 

the algebra still being of type VI/j. 

The further rescaling (12.21) and the limit D ^ transform (12.24) into a Bianchi 
type V algebra, but /c^^^ has to be replaced by 

^(3) = (^r3) + 1^(1)) ------> {n/m)6t + {H/m)e-^''''5t. (12.25) 

When J = on top of ^ 0, the Bianchi type reduces to I. 

The subcase with g > j^/4 (Bianchi type Yllh with the free parameter j/ {Ag — j'^Y^'^) 
is given by eqs. (10.18) - (10.19) in Ref. 3, with yo = 0. It is transformed back to the 
Plebahski coordinates by the inverse of (10.8) there, and the result is very similar to our 
(12.10) - (12.11). Only the definitions of U, V and y, and a few signs in the metric are 
different: 

U ■=Mcos{DY) + iVsin(DF), V := M sm{DY) -Ncos{DY), 
P = e-^-^/^C/, y = -^P + -e-^^/V, 

2^ g 

911 = 2/' + h^iP^ + 2Pe-^'^t/y/ii3 + (7r/i^)e-^'^/i33, 
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gi2 = gthii + /ii2 + Vgt{V/U + -j/D)hn + {rV/U)h23 

JV 



2DU 

913 = hsP + re-^'^/V/i33, 
922 = + 29thi2 + {jT/D)gt{9this + h^s) + /i22 + 9{^9t/Dfhss 

923 ^ P'^ {gthi2, + /l23 + ^^^33) , 5-33 = ^33- 

The hij are here functions of the argument 

T = t z. 

^9 



;i2.26) 



;i2.27) 



The redefinitions in the constants and functions needed here are again (12.16) - (12.17) 
together with 

/ii2 = -gTHn/cuo'' + Hi^/coo - |^^i3 - r'^r/i33, 
h22 = -{9TfHn/u;o^ - 2gThu + H22 - {jVg/D) (gT^H^^/uj^ + ^23) - g{TgT/Dfh^^, 



h 



23 



-gTHrs/uJo + H23 ^tt^'^ss- 



2D 



;i2.28) 



The metric resulting after all the redefinitions is 



L>(m2 + n2) ' 

511 = i^oyf + HnP^ + Te-^^ [2UVH,s + {h/D)^^ 
gu = P [(D/r)zHn + + (j/2 + DV/U)zH,s 

+ {rV/U)H23 + t(D + \jV/U^ zhs3^ , 

^13 = //l3P + fe-^'^/'y/l33, 

g22 = {Dz/tfHu + iD/t){2zHn+jz''Hi3) + H22+3ZH23 + gz'^h^z, 

g23 = {D/f)zH^^ + H23 + (j72)2:/i33, ^?33 = h^3 = ^33, (12.29) 
Just as before, in the limit ujq ^ Q the Hij will depend only on t, and the shearfree 



hmit is found in the same way: Hij — —CijR , C33 = 1. 

The k — —1 Priedmann limit is now obtained in two ways: either 

h = HD, D^Q, 

or 

Cii = Ci3 = 0. 



;i2.30) 



(12.31) 
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In the first case the Priedmann hmit is 



-R" [e-^'^/'(mCi3 + nV)dx + (C23 + 3z/2)dY + d^l ' , 



(12.32) 



where 




(12.33) 



ds^ = dt^ - R\t) [D^e-^^/^Udx + (D12 + Dz)dY] - {D22RdYf 
-R^ \e-^^l'^VVdx + (C23 + jz/2)dY + d^l^ 



(12.34) 



where Dn is defined as in (12.33), but with Cu = C13 = 0. 

The KiUing fields before redefinitions are still given by (12.13), but, in consequence of 
the different definitions of P and y in the present subcase, the Bianchi type is VII/j. In 
the Friedmann limit defined by (12.31), the Killing fields are transformed only by (12.16) 
- (12.17) followed by ujq 0, and they still form a type Vllh algebra. When (12.30) is 
imposed on top of (12.16) - (12.17) and uq 0, the Killing fields become the same as the 
limit D — > of (12.24) - (12.25), i.e. the Bianchi type becomes V. This is an illustration 
of the fact, mentioned in sec. I, that the k — —1 Robertson- Walker geometry is a subcase 
of two Bianchi types simultaneously, they are exactly V and Vllh- 

In both cases, the A; = Friedmann limit follows from (12.32) and (12.34) when j = 0. 
In the first case, the algebra of the KiUing fields becomes type I, in the second case it 
becomes type VIIq, which is another illustration of the same kind of duality. 

Finally, the third subcase of case 2.2.1.2 is given by eqs. (10.20) - (10.21) in Ref. 3, 
with Ho — 0. There is one more typo there, the correct formula for y is 



This one is of Bianchi type IV. When transformed back to the Plebahski coordinates (by 
the inverse of (10.8) in Ref. 3), it becomes 



y = -2P/j - AMe-^^/'^/f. 



(12.35) 



p = e-^^/'^{MY + N), 
gii = + hnP^ - 2(7/M)e-^'^/'P/ii3 + /Mfe'^^'h^^, 



8M2 



fit 



h: 



'33 



) 
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/ 1 j^jt \ 

923 = ( -^fthis + h23 + ■^j^h^sj ' 



933 — ^33- 



where the hij are arbitrary functions of 



t- 



;i2.36) 



;i2.37) 



The redefinitions needed to calculate the limit ujq ^ are 

(M, N, 7) = (m, n, h)uJo, hu = Hu/ujq^, hi3 = Hi^/uq, 



h 



1 j^/iT 



h22 = -(j^r/4)2i/n/^o' - ^fTh 



12 



His + 



22 



h 



23 



23 



fhT 

-]fTHn/uJo + H23- 



^33, 
fhT 



k 



33- 



(12.38) 



We will denote, as before, {y,P) = WQ{y,P), and choose Y as the new x^-coordinate, 
so that dy = cuoPdY. The metric that results is 

911 = i^oyf + HnP^ - 2{h/m)e-'^/^PHn + (h/mfe-^^ hss, 



gu = P {{myh)zHn + + {j/2)zHis 



-(mY + n)~'^ 



mzHis + {h/m)H23 + ^{jh/m)zh33 



}• 



(7l3 = i/l3P-(V^)e-^'^/'/l33, 

922 = {m'z/hfHii + 2{myh)zHi2 + {jrif/h)z^His + H22 + 3ZH23 + {jz/2fh33, 

923 = {myh)zHi3 + H23 + {j/2)zh3s, ^33 = ^^33 = ^33, (12.39) 



In the limit a;o — * 0, the Hij will depend only on t. 
The A; = — 1 Friedmann limit is now obtained when 

Hij^ -CijR^it), h^Hm, 

The k — limit will result when j = in addition. 



m — > 0. 



;i2.40) 
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XIII. The cases 2.2.2 of Ref. 3. 

The case 2.2.2.1.1 is given by eqs. (11.11) - (11.12) in Ref. 3. The transformation back to 
the Plebanski coordinates is the inverse of (11.10), and the transformed metric is 

gu = y'^hu, gi2 = {jt/o){l - /ill) + hi2 - Ajthis, 

9i3 ^ B{j + a)yhis, 



922 = ( —] {hu - 1) - 2^/ii2 + 2A^his + h22/y^ - 2{Ajt/y^)h23 + {Ajt/yfhss, 
ay I ay^ ay^ 



(jt Ajt \ 
hi3 + h2z/y /i33 , Qzz^ B'^{j + afh^^, (13.1) 
ay y J 

where a, A, B and j are arbitrary constants, and the hij are arbitrary functions of the 
argument 

T^t + ^^^z. (13.2) 
Aa ^ ' 

The KiUing fields for the metric (13.1) are 

^(1) ~ ^(2) ~ ~ y^2i 



i^{3) — y 



B{j/a + 1)6^ - Bj{ay)-^S^ - AS^ , (13.3) 



and they form a Bianchi type VI/i algebra, with the free parameter {j — a)/ {j + a) . The 
analog of (5.1) is 

k^^ - Bj{ay)-^kl^ = y-^l'^ [B{j/a + 1)^" - {A/n)w^] , (13.4) 
and the King-Ellis measure of tilt is 

^/^u"N^ = Ay^-^'". (13.5) 
By a simple transformation of the 2;-coordinate we can achieve the same result as if 

S(j + a) = l, (13.6) 

and we will assume this now. 

The redefinitions needed to calculate the limit tuo — >^ are 

y = uoy, A = A/uq, hu = Hu/uq^, hi3 = His/ujq, 

hi2 = -{j/a)T + {jT/a)Hn/uo^ + H^/uJo + AjTHi^/uo'', 
h22 = (jT/af (l - i^ii/a;o') + 2{jT/a)hi2 - 2{AfTya)H,3/uJo' + H22 
+2AjTh23/u;o - {AjTfhss/coo'', 

h23 = {jT/a)H^s/uJo + H23 + AjThs/cuo. (13.7) 
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The redefined metric is 



9n 



922 



9l2 = 



a?A 

2 



z{Hii - loq^) + Hi2 + {j/a)zHi3, gis = yHis, 



+2{j/a)zH23 + Uz/afh33 



2^^^^i3 + H22 



923 = y 



J 



■zHi3 + H23 + {j/a)zh3s 



933 — ^33 — -f^33) 



.a A J 
In the hmit a;o ^ 0, the Hij will depend only on t. The A; = — 1 Priedmann limit will then 
result when 

Hij = -CijR (t), j = -a, 
The /c = Priedmann limit will result when 



00. 



;i3.9) 



Hij — —CijR^if), 



Cl3 = J = 0, 
e^", A; ^ 0. 



C22 = {D22/kf, 



y = e'^", k^[). (13.10) 

The limit (13.9) transforms the Killing fields (13.3) into a Bianchi type V algebra, 
provided k'^r^-^ is redefined to k'^^ = —{ujQ^^"'/A)k'^r^^, and the limits A ^ 00 and j = —a are 
tuned so that A{j + a) — 00 (for example, A = a/{j + a) and a — >• 00). The limit (13.10) 
will transform (13.3) into a type I algebra, but fc^-^ has to be redefined as above, and in 
addition k^2) has to be redefined to k'^^ — kk^2)- 

The formulae for the case 2.2.2.1.2 simply follow from those above. This case has a 
completely different outlook only in the coordinates adapted to the Killing fields that were 
used in Ref. 3. When transformed to the Plebanski coordinates, it becomes the subcase 
of (13.1) given by 

B{j + a) = 1, J = -a, A = Ai/a, (13.11) 

where the Ai defined above stands in place of the A from eqs. (11-17) in Ref. 3. Then the 
redefinitions needed are (13.7) with Ai — A/cuo and j — —a, and the redefined metric is 

gn = fHn, 9i2 = -{z/A){Hn - uJo^) + H^ - zHi^, g^s = yHi^, 



922 = y 



-.-2 



{z/A)\Hii - u;o') - '2{z/A)Hi2 + 2{z''/A)H^3 + H22 - 2^/^23 + z'h. 
923 = y~\-{z/a)Hi2, + H23 - zhz3], 933 = ^33 = -f^33, 



(13.12) 

The k — —1 Priedmann limit results now by (13.9) when ujq = 0, and the A; = Priedmann 
limit results from (13.12) when 

ujo = 0, H,j = -CijR\t), {y, z) = (e'=^ e^^), 

{H^2, i/13) = (^12, Gis)/k, {H22, H23, i/33) = (^22, ^23, Gss)/e, k ^ 0. (13.13) 
The KiUing field k^^-^ is different here 



k 



(3) 



1/5% -Ini/^i- ^1/^3, 
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;i3.i4) 



while the two others are as in (13.3). The case 2.2.2.1.2 required a separate consideration 
in Ref. 3 only because of the logarithm term in the Killing field. In calculating the 
limits a;o — >■ and A — >■ oo, this vector field has to be redefined similarly as before. For 
calculating the A; = Priedmann limit, ^"3^ and k'^^} ^^^e to be redefined by {k'^-^, k'^^-^) — 
k{k^,),kf,^) ________> {-S-„-AS-s). 

Finally, the case 2.2.2.2 from Ref. 3 (eqs. (11.18) - (11.27)) is of Bianchi type I, with 
the velocity field being tangent to the symmetry orbits, so it has no Friedmann limit at 
all. 

XIV. Summary. 

All the metrics derived in Refs. 1-3, that correspond to rotating hypersurface-homogeneous 
dust models, have been checked here for the existence of a Friedmann limit. It was found 
that such a limit exists for all those cases listed in Refs. 2 and 3, where the matter-density 
is not constant along the fiow. However, in at least one class (see sec. Ill), the Friedmann 
model will have no rotating parent solution, but will instead be a separate subclass. 

Along the way, the nonstationary metrics were all transformed to such a form, in which 
the limit of zero rotation can be explicitly calculated. The transformation/reparametrization 
leading to this form is nonsingular and invertible in each case, but it becomes singular when 
uj 0. The limits c<j = all have nonzero shear. Thus, a whole collection of metrics gener- 
alizing those of Friedmann was found that can be used in studying spatially homogeneous 
exact perturbations of the latter. 

The Class A Bianchi-type metrics (those in which the structure constants have the 
property C"ac = 0) are known to admit a lagrangian/hamiltonian formulation^^. Those of 
them that obey the Einstein equations with a rotating dust source (types VIq, VIIq, VIII 
and IX) were studied by Ozsvath^'^^. The lagrangians and hamiltonians were explicitly 
found in Refs. 7 and 13, and the Einstein equations in the hamiltonian form were then 
transformed to such variables, in which they become analytic. This should prove the 
existence of solutions. 

Two more papers, specifically devoted to rotating spatially homogeneous dust solutions, 
are those of Bchr^^ (where a subclass of type IX models was investigated) and of this 
author"^ (discussing a subclass of type V models). In both of these, the Einstein equations 
were transformed, simplified, investigated for known limiting cases and for Lie symmetries, 
but no explicit solutions were found. A (hopefully) complete overview of other solutions 
with rotating matter source is given at the end of Ref. 3. 

It is hoped that the present paper will be helpful in picking out those models for future 
investigation that promise interesting physics or geometry. 
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CAPTION TO THE DIAGRAM 

The diagram shows how the different Bianchi types can be specialized by taking the 
zero limit of one or more of the structure constants. This allows to recognize (by the rules 
given at the end of sec. I), which Friedmann models can possibly be contained as limits in a 
given class - see text. All the possibilities are actually realized in the collection considered 
in the paper. 
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